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Considerable effort has been applied towards understanding the precise shaping mechanisms responsible for the diverse range of 
morphologies exhibited by planetary nebulae (PNe). A binary companion is increasingly gaining support as a dominant shaping 
mechanism, however morphological studies of the few PNe that we know for certain were shaped by binary evolution are scarce or 
biased. Newly discovered binary central stars (CSPN) from the OGLE-III photometric variability survey have significantly increased 
the sample of post common-envelope (CE) nebulae available for morphological analysis. We present Gemini South narrow-band 
images for most of the new sample to complement existing data in a qualitative morphological study of 30 post-CE nebulae. Nearly 
30% of nebulae have canonical bipolar morphologies, however this rises to 60% once inclination effects are incorporated with the aid 
of geometric models. This is the strongest observational evidence yet linking CE evolution to bipolar morphologies. A higher than 
average proportion of the sample shows low-ionisation knots, filaments or jets suggestive of a binary origin. These features are also 
common around emission-line nuclei which may be explained by speculative binary formation scenarios for H-deficient CSPN. 
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£ — , ■ 1. Introduction thought to produce bipolar nebulae and a mild one elliptical neb- 

O ' uiae ( Livio & Soker 1988; Bond & Livio 1990). The complexity 

^ ■ A single shaping mechanism behind the extraordinary morpho- of thefje models Hmits their predictive power when it comes t0 

, logical diversity of planetary nebulae (PNe) remains elusive even morphological studies of observed nebulae, however there are 

• , after more than 20 years of discussion (Bahck & Frank 2002). SQme indications that the density contrast could be as high as 

. £ ; Highly axisymmetric, non-spherical forms in particular are very 5;1 (Sandquist et aL 1998) . This remains t0 be proven with more 

V , difficult to explain without a binary engine that is capable of sup- advanced simu i ations , but if such a density contrast is accurate, 

g • plying sufficient angular momentum to heavily shape the neb- then we may suspect tQ find post . CE nebulae dominated by bipo . 

CZ | ula (Nordhaus, Blackman & Frank 2007). The binary hypoth- ]ar morp hologies 

esis encapsulated by De Marco (2009) simply states that PNe 

derive from binary progenitors more easily than from single pro- It may therefore be surprising that Bond & Livio (1990) 
genitors. Gathering the direct evidence of binary central stars did not find a large proportion of bipolar nebulae amongst -12 
(CSPN) required to test the binary hypothesis has however been post-CE nebulae wherein only NGC 2346 satisfied the strict 
a slow process with only around a dozen known until recently morphological criteria of canonical bipolar nebulae (Corradi & 
(De Marco, Hillwig & Smith 2008) compared to -2600 PNe Schwarz 1995). With such a small sample it is difficult to de- 
known in the Galaxy (Parker et al. 2006; Miszalski et al. 2008a). termine whether the theoretical expectations or the morphologi- 
Such a small sample has considerably hampered attempts to cal classifications are responsible. We suspect projection effects 
reveal morphological trends amongst the nebulae surrounding are primarily responsible for the apparently low bipolar frac- 
binary CSPN. We know their nebulae must have been ejected tion amongst post-CE PNe. Few morphological studies have at- 
during the common-envelope (CE) phase (Iben & Livio 1993) tempted to account for the complex 3D structure of each nebula 
and attempts to model this process have suggested the geom- when projected onto the 2D plane of the sky (e.g. Manchado et 
etry of the mass loss is subject to a density contrast between al. 2004). The appearance of each nebula is therefore strongly 
the equatorial and polar directions (Livio & Soker 1988; Bond dependent on its inclination angle with respect to the plane of 
& Livio 1990; Sandquist et al. 1998). A high density contrast is the sky (e.g. Frank et al. 1993) and recovering the instrinsic 3D 
morphology requires detailed kinematic modelling. Consistent 



* Based on observations made with Gemini South under programs application of these techniques to post-CE nebulae has only just 
GS-2008B-Q-65 and GS-2009A-Q-35. commenced with Mitchell et al. (2007) revealing an intrinsic 
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bipolar morphology to A 63. Mitchell et al. (2007) demonstrated 
effectively the utility of knowing the orbital inclination of the bi- 
nary CSPN in building a comprehensive model of the nebula in 
conjunction with kinematic data. 

Even without kinematic information there is still room for 
qualitative classification within the framework of geometric 
bipolar models forpost-CE nebulae. Pollacco & Bell (1997) con- 
structed such a model and made favourable comparisons with A 
41, A 46, A 63 and A 65 suggesting they are bipolar nebulae. 
Similarly, De Marco (2009) considered inclination effects with- 
out the assistance of models to estimate as much as 70% of post- 
CE nebulae could indeed be bipolar. This estimate is however 
disputable in that it was heavily biased by a small sample and in 
some cases the highest quality images available in the literature 
were not used (e.g. A 65). A detailed reappraisal of a large sam- 
ple of post-CE nebulae is therefore required to substantiate the 
suspected high bipolar fractions from recent observations and to 
compare them with current theoretical expectations. 

To conduct this reappraisal we combine the small sample of 
prevously known post-CE nebulae (De Marco et al. 2008) with 
those found by Miszalski et al. (2008b; 2009, hereafter 'Paper 
F). In Paper I we cross-referenced previously known PNe and 
those from the Macquarie/AAO/Strasbourg Ha (MASH) PNe 
catalogues (Parker et al. 2006; Miszalski et al. 2008a) with high 
quality time-series photometry from the OGLE-III survey to- 
wards the Galactic bulge (Udalski 2009). Our comprehensive 
approach yielded a more representative sample of 21 new peri- 
odic binary CSPN ideally suited to a less-biased morphological 
study. 

This paper is structured as follows. Section [2] discusses our 
sample selection and available observations. Geometric models 
and a classification scheme are introduced in Sec. [3] These are 
used to guide morphological classification of the sample in Sec. 
|4] where individual PNe are discussed in detail. Section [5J] de- 
scribes the frequency of the morphological types, while sections 
15.21 and 15.31 discuss low-ionisation structures. We conclude in 
Sec.lH 



2. Observations 

2. 1 . Sample selection 

We are solely concerned with the spatially resolved nebulae sur- 
rounding periodic close binary central stars. Spectroscopic con- 
firmation of each central star is required to construct a clean 
sample, however those without confirmation may also be accom- 
modated depending upon the quality of their candidate CSPN. 
Central stars tabulated by De Marco (2009) that appear to be 
visual binaries, have composite spectra or main-sequence near- 
infrared colours are omitted from our sample since their periods 
are unknown. Visual binaries are not directly comparable to our 
sample because their influence on nebular shaping is expected to 
be less dramatic than in close binaries (Soker 1997), whereas un- 
confirmed composite or cool central stars without periods may 
just be superpositions (Mendez 1989). 

Table [T] lists 15 nebulae from De Marco et al. (2008) and 
18 nebulae from Paper I that meet our selection criteria. These 
samples will be referred to as the known sample and the OGLE 
sample, respectively. Column three gives the orbital period, col- 
umn four the major and minor axis diameters, and column five 
indicates whether stellar eclipses occur in the lightcurve. Other 
columns are reserved for Spitzer IRAC photometry (Sec. 12. 3i or 
morphological classification purposes (Sec. 13. It . 



The binary status of 14 members from the known sample 
is well established with plenty of observational evidence (De 
Marco et al. 2008). Only periodic variability is reported for HaTr 
4 (Bond & Livio 1990), but we include it in our sample based on 
the well-centred position of the CSPN. We exclude SuWt 2 from 
De Marco et al. (2008) because the inclination of the bipolar 
nebula deduced from the equatorial ring is inconsistent with the 
double A-type eclipsing central star and the true central star may 
yet turn out to be one of two other candidates revealed in high 
quality imaging (Smith et al. 2007). Also omitted is Hb 12 whose 
period is purportedly irregular (Hillwig, priv. communication to 
De Marco 2009), as well as the group of A 35-type binaries LoTr 
1, LoTr 5 and A 35 since the PN status of A 35 is disputed (Frew 
2008; De Marco 2009). 

We conducted spectroscopic observations of 14 binary 
CSPN candidates from Paper I under Gemini programs GS- 
2008B-Q-65 and GS-2009A-Q-35. The Gemini Multi-Object 
Spectrograph (GMOS, Hook et al. 2004) was used in its longslit 
configuration mostly with the B1200 grating to target the pri- 
mary component over the wavelength range 4085-5550 A. 
Future papers will discuss the spectra in detail, however we can 
discuss here those either not confirmed, awaiting confirmation or 
having some other significance. Most candidates have been con- 
firmed after identification of features belonging to the primary 
(HI and Hell absorption) or the irradiation of the secondary (CIII 
and Mil stellar emission, e.g. Pollacco & Bell 1993, 1994). For 
example, features of both components were detected in K 6-34, 
leading us to revise downward its period to 0.20 days. No fea- 
tures of the primary were visible for PHR 1744-3355 and PHR 
1801-2718, each of which shows G-type features shifted by at 
least 50 km/s from the nebula. Although we cannot rule out bi- 
nary motion as the cause of the radial velocity difference, both 
associations seem unlikely at this stage and so we exclude them 
from our sample. A composite spectrum is however detected in 
PHR 1804-2913, but radial velocity monitoring will be required 
to confirm whether its unusual 6.7 day period is indeed an orbital 
period. Further additions to our sample from Paper I include H 
1-33, Bl 3-15 and PHR 1804-2645 that have well-located binary 
CSPN candidates that were too faint or too reddened to show any 
stellar features in their AAOmega spectra (described in Paper 
I). Another member of this group would be JaSt 66, but lack- 
ing any resolved imagery we have omitted it from our sample. 
The remainder of our sample includes PHR 1759-2915, MPA 
1759-3007 and PHR 1801-2947 whose binary CSPN candi- 
dates are yet to be confirmed, however these will be treated sep- 
arately from the other objects given their higher uncertainty. 

2.2. Narrow-band optical imaging 

As part of the acquisition process during our Gemini South pro- 
grams we took advantage of the exceptional imaging capabilities 
provided by GMOS (Hook et al. 2004). Table |f] records the fil- 
ters and exposure times of images used in this work. The central 
CCD was read out with 2x2 binning to give a 5.5 x 2.5 arcmin 
field-of-view sampled at 0.145"/pixel. The central wavelength 
and FWHM of the filters are 499.0/4.5 nm (OIII), 514.0/8.8 nm 
(OIIIC), 656.0/7.2 nm (Ha), 662.0/7.1 nm (HaC) and 672.0/4.4 
nm (SII). The Ha filter includes Ha and both [Nil] lines, while 
the HaC filter effectively acts as an [Nil] /I6584 filter. 

Figure Q] presents a montage of mostly Ha+ [Nil] images 
of the OGLE sample (Tab. Q]). The main contribution comes 
from our Gemini programs, followed by the ESO NTT/EMMI 
programs 079.D-0764(B) and 67.D-0527(A) for H 1-33 and 
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Table 1. PNe containing close binary central stars whose nebular morphologies are assessed in this work. 
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Period 


Diameter 


Eclipsing 
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Model// 


Morphology 


Traits 


IRAC 


Ref. 






(days) 


(arcsec) 




(deg) 








12 3 4 




005.0+03.0 


Pe 1-9 


0.14 


12 x 14 


Y 


-90 


_ 


MS 


L 


NUNU 


b 


355.3-03.2 


PPA 1747-3435 


0.22 


12 


N 


_ 


_ 


Irr 


_ 


UUUU 


a 


355.7-03.0 


H 1-33 


1.13 


5 x 12 


N 


_ 


A/45 


CB 


G 


U YU Y 


a 


354.5-03.9 


Sab 41 


0.30 


95 x 194 


N 


_ 


A/40 


LB 


L,J,R 


UUUU 


b 


000.6-01.3 


Bl 3-15 


0.27 


4x11 


N 


_ 


A 


LB 


G,J 


Y Y Y Y 


a 


359.1-02.3 


M3-16 


0.57 


9 x 13 


Y 


-90 


_ 


Irr 


L,J 


Y Y YY 


b 


357.6-03.3 


H2-29 


0.24 


11 x 14 


Y 


-90 


A/90 


CB 


L,R 


UUUU 


b 


000.2-01.9 


M2-19 


0.67 


7 x 19 


N 


_ 


A/45 


CB 


G,R 


Y Y YY 


b 


358.7-03.0 


K6-34 


0.20 


14 


N 


_ 


_ 


MS/LB 


G,L,J?,R 


NNN Y 


b 


357.0-04.4 


PHR 1756-3342 


0.26 


22 


N 


- 


B/30 


PB 


- 


UUUU 


b 


001.8-02.0 


PHR 1757-2824 


0.80 


4 x 18 


Y 


-90 


A/90 


CB 




NNNN 


b 


001.2-02.6 


PHR 1759-2915 


1.10 


25 x 30 


Y 


-90 


A/90 


CB 


_ 


NNNN 


b 


000.5-03.1a 


MPA 1759-3007 


0.50 


17 x42 


N 


- 


A/55 


LB 


- 


UNUN 


b 


001.9-02.5 


PPA 1759-2834 


0.31 


17 x 35 


N 




B/45 


LB 


I 


NNNN 


b 


357.1-05.3 


BMP 1800-3408 


0.14 


30x55 


Y 


-90 


- 


Irr 


I 


UUUU 


b 


000.9-03.3 


PHR 1801-2947 


0.32 


34 


N 




B/30 


PB 


_ 


NNNN 


b 


001.8-03.7 


PHR 1804-2913 


6.66 


3.5,8 


N 


_ 


_ 


MS 


I 


UNNN 


b 


004.0-02.6 


PHR 1804-2645 


0.62 


24x45 


Y 


-90 


A/90 


CB 


G 


NNN Y 


b 


144.8+65.8 


BEUMa 


2.29 


180 


Y 


84 


A/85 


LB 


_ 


UUUU 


J> n 


055.4+16.0 


A 46 


0.47 


53 x 70 


Y 


80 


A/80 


LB 


I 


UUUU 


b,c,g 


053.8-03.0 


A 63 


0.46 


42 x 290 


Y 


88 


A/90 


CB 


J 


UUUU 


g 


136.3+05.5 


HFG 1 


0.58 


540 


N 


28 


_ 


Irr 


LL 


UUUU 


f 


253.5+10.7 


K 1-2 


0.68 


50 x 130 


N 


-45 


A/45 


LB 


L.J 


UUUU 


b,i,m 


283.9+09.7 


DS 1 


0.36 


290 x 335 


N 


62 


_ 


PB 


LL 


UUUU 


b 


329.0+01.9 


Sp 1 


2.91 


120 


N 


_ 


A/5 


LB 


G,R 


Y Y Y Y 


c 


349.3-01.1 


NGC 6337 


0.17 


60x100 


N 


9 


A/10 


LB 


G,L,J,R 


Y Y Y Y 


1 


005.1-08.9 


Hf 2-2 


0.40 


25 


N 






MS 


L,R 


UUUU 


b.d 


017.3-21.9 


A 65 


1.00 


220 x 240 


N 




A/10 


LB 


J 


UUUU 


c.h 


009.6+10.5 


A41 


0.23 


15 x 22 


N 


66 


A/70 


CB 


R 


UUUU 


b.g 


335.2-03.6 


HaTr4 


1.74 


24 x 30 


N 




B/45 


PB 




UUUU 


c 


215.6+03.6 


NGC 2346 


15.99 


58 x 172 


N 


50 


A/50 


CB 


G,R 


Y Y Y Y 


c 


135.9+55.9 


SBS 1150+599A 


0.16 


4 


N 






Irr 


L? 


Y Y Y Y 


k 


341.6+13.7 


NGC 6026 


0.53 


60 


N 


82 




Irr 




UUUU 


d 



Morphologies: CB: Canonical bipolar; LB: Likely bipolar; PB: Possible bipolar; MS: Multiple shells; Irr: Irregular 
Morphological Traits: G: Gatley's rule satisfied; I: ISM interaction; L: Low-ionisation knots or filaments; J: Jets; R: Equatorial ring 
IRAC: U: Unavailable; Y: Nebular detection; N: Nebular non-detection. 

References: (a) Paper I, (b) this work, (c) Bond & Livio (1990), (d) Schwarz, Corradi & Melnick (1992), (e) Walsh & Walton (1996), 

(f) Heckathorn, Fesen & Gull (1982), (g) Pollacco & Bell (1997), (h) Hua, Dopita & Martinis (1998), (i) Corradi et al. (1999), (j) Bond (2000), 

(k) De Marco (2009), (1) Corradi et al. (2000), (m) Exter, Pollacco & Bell (2003), (n) Liebert et al. (1995) 



Table 2. GMOS acquisition images used in this work. 



Name 


Filter 


Exp. 


Name 




Filter 


Exp. 


Pe 1-9 


Ha 


180 


K6-34 




Ha 


240 




OIII 


120 






OIII 


240 




OIIIC 


120 






OIIIC 


120 


Sab 41 


OIII 


240 


PHR 1756- 


-3342 


Ha 


300 




Ha 


240 






OIIIC 


120 




HaC 


120 


PHR 1757- 


-2824 


Ha 


150 


M3-16 


Ha 


120 






HaC 


60 




OIII 


120 


PPA 1759- 


2834 


OIII 


300 


H2-29 


Ha 


60 






g' 


90 




SII 


60 


BMP 1800 


-3408 


Ha 


300 




OIII 


60 






r' 


90 


M2-19 


OIII 


120 


PHR 1804- 


-2913 


Ha 


60 




Ha 


120 






HaC 


60 




SII 


120 











Bl 3-15, respectively (see Paper I), and the AAO/UKST 
SuperCOSMOS Ha survey (SHS, Parker et al. 2005). We fur- 
ther add Spitzer colour-composite images described in Sec. 



12.31 Image dimensions are all 30 x 30 arcsec 2 except for 
PHR 1804-2913 (15 x 15 arcsec 2 ), PHR 1804-2645, MPA 
1759-3007, BMP 1800-3408, PHR 1801-2947 and PHR 
1759-2915 (1 x 1 arcmin 2 ), and the SHS quotient image of 
Sab 41 (4 x 4 arcmin 2 ). Throughout this work all figures are 
presented in the standard orientation unless otherwise indicated 
(North is up and East is to left) and each image is centred on 
the central star. The SHS images are also complemented later 
by broad-band B j images from the SuperCOSMOS Sky Surveys 
(SSS, Hambly et al. 2001) in constructing colour-composite im- 
ages from SHS Ha (red), SHS Short-Red (green) and SSS B, 
(blue). 



Our primary focus in this work is on the new OGLE sam- 
ple rather than previously well-studied known sample. Section 
|4]will present a selection of more detailed images of the OGLE 
sample than that provided in Fig. [TJ A reappraisal of available 
images for certain members of the known sample will also be 
made alongside the OGLE sample, however for those not cov- 
ered we refer the reader to the final column of Tab.Q] 
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Fig. 1. A montage of mostly Hci'+[NII] false-colour images and 
Spitzer IRAC colour-composite images of the 18 nebulae se- 
lected from Paper I for morphological analysis. 



2.3. Gatley's rule and mid-infrared imaging 

As a complement to narrow-band optical imagery, near- and 
mid-infrared imagery can help reveal an intrinsic bipolar mor- 
phology regardless of the projected 2D morphology. This devel- 
opment started with Zuckerman & Gatley (1988) and was later 
solidified by Kastner et al. (1994, 1996) who conducted narrow- 
band near-infrared (NIR) imaging surveys of PNe in the 2.122 
pm H2 emission line. They found a strong propensity for canon- 
ical bipolar PNe to display H2 emission in their waist regions. 
Excluding canonical bipolar PNe, H2 detections were found only 
amongst torus or ring-like structures interpreted to be inclined 
equatorial waists of bipolar nebulae. This is known as Gatley's 
rule, namely that the detection of H2 emission is a unique marker 
of a canonical bipolar PN (Kastner et al. 1996). 

While we acknowledge that Gatley's rule has yet to be rig- 
orously tested with large and unbiased samples, there are some 
consistent observational and theoretical indications that the rule 
holds and is suitable for use in identifying bipolar PNe. An inde- 
pendent study by Guerrero et al. (2000) found definite H2 de- 
tections in 13/15 bipolar PNe. These encouraging results are 
supported theoretically by Aleman & Gruenwald (2004) who 
modelled the physical conditions responsible for the presence 



of H2 in PNe. They found higher H2 mass fractions are associ- 
ated with higher central star temperatures, which is consistent 
with the higher than average temperatures amongst bipolar PNe 
(Corradi & Schwarz 1995), thereby providing a simple physical 
explanation of Gatley's rule. 

Although we lack narrow-band near-infrared H2 imaging for 
our sample, we can make use of the next best thing, observations 
made with the Infrared Array Camera (IRAC; Fazio et al. 2004) 
onboard the Spitzer Space Telescope. The four IRAC bands are 
centred on 3.6 pm (band 1), 4.5 pm (band 2), 5.8 pm (band 
3) and 8.0 pm (band 4). Band 4 in particular serves as a good 
tracer of H2 emission containing two H2 lines and an [Ar III] 
line (Hora et al. 2004). Bands 2 and 3 are also sensitive to H2 
emission but not as much as band 4. It follows that a strong band 
4 detection and moderate to absent detections in the other bands 
satisfies Gatley's rule, i.e. the PN is highly likely to have an in- 
trinisc bipolar morphology. Emission lines from atomic transi- 
tions and PAHs do occur within the broad IRAC bands however 
(e.g. Cohen et al. 2007), and while their application in this con- 
text is not free from ambiguity, it nonetheless serves as a useful 
tool in the absence of proper narrow-band H2 imaging. Indeed, 
Kwok et al. (2008) found excellent agreement between a narrow- 
band H2 image of NGC 6583 and the Spitzer IRAC band 4 image 
while demonstrating the bipolar (double-coned) morphology of 
NGC 6583. Another dramatic case where H2 emission unam- 
biguously identifies bipolar morphology not necessarily evident 
in optical emission lines can be seen in the velocity integrated 
H 2 image of NGC 7027 (Cox et al. 2002). 

We have searched the Spitzer science archive using the leop- 
ard tool for publicly available IRAC observations for our sam- 
ple. Tab. Q] records the status of each IRAC band as either be- 
ing unavailable (U), a nebular detection (Y) or a non-detection 
(N). Almost all observations were conducted as part of the ongo- 
ing GLIMPSE surveys (Benjamin et al. 2003; Churchwell et al. 
2009). Since GLIMPSE is confined mostly to \b\ < 1°, much of 
the known sample at higher latitudes have no available images. 
Images of SBS 1150+599A were taken serendipitously during 
observations of the blue compact dwarf galaxy SBS 1150+599 
(Wu et al. 2006), whereas NGC 2346 was observed as part of 
program ID 68 (PI: Fazio). All colour-composite Spitzer images 
presented in Fig.Q]and elsewhere in this paper were constructed 
using IRAC bands 4 (red), 3 (green) and 2 (blue). In these im- 
ages a ruddish or red-orange hue, expected to coincide with the 
waist of a bipolar PN, satisifies Gatley's rule. Non-detections for 
the faintest MASH PNe could be explained by their very low 
surface brightness. 

3. Classification Scheme 

3. 1 . Geometric models 

We have constructed two rudimentary geometric models for a 
range of inclination angles to compare qualitatively with the pro- 
jected 2D morphologies of our sample. Columns seven and eight 
of Tab. [1] encapsulate these comparisons in the rough choice of 
model and inclination, and the overall morphology, respectively. 
Depending on the inclination of the system and which geometric 
model is used, the overall morphology may be one of canoni- 
cal bipolar (CB), likely bipolar (LB) or possible bipolar (PB). 
Nebulae showing multiple-shell (MS) or irregular (Irr) mor- 
phologies are classified separately without the models. The mor- 
phological traits listed in column eight of Tab. Q] show whether 
Gatley's rule is satisfied (G, see Sec. \2.3i . ISM interaction is 
likely (I), or an equatorial ring is present (R). Low-ionisation 
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structures (L) and jets (J) are defined in Sec. 13.21 and discussed 
further in Sec. 15.21 

Figure|2]presents our slight extension of the Pollacco & Bell 
(1997) model (hereafter, Model A) which has been used success- 
fully to reveal the intrinsic bipolar morphology of A 63 (Mitchell 
et al. 2007). The components are a main outflow (blue hyper- 
boloid), an equatorial ring of low-ionisation material (large red 
torus), polar ejecta or jets (small red torii) and a central star 
(black dot). Note at low inclination angles the projected radial 
distance of the jets depends heavily on their height above the 
equatorial plane. We do not demand a precise match to Model 
A since the aspect ratios or even overall shape of the main out- 
flow are considered to be quite flexible. Other shapes may be 
incorporated including a 'peanut' or double cone (Kwok et al. 
2008) structure, but these are less common amongst the ob- 
served morphologies. Matches to Model A are typically given 
the likely bipolar (LB) classification or alternatively canonical 
bipolar (CB) if they meet the strict criteria outlined in Sec. 14. II 

Also included in Fig. [2]is a more rudimentary version of the 
Frank et al. (1993) model (hereafter, 'Model B'). Model B con- 
sists of a main outflow (blue modified prolate spheroid) and a 
central star (black dot), but it could conceivably have the same 
low-ionisation ring or jets as seen in Model A. The formation 
of Model B involves a much lower density contrast between the 
equatorial and polar directions than Model A. It would there- 
fore be considered more elliptical than bipolar, especially since it 
lacks the conventional 'pinched' or concave waist. This conven- 
tion however does not account for evolutionary and inclination 
effects befitting bipolar nebulae (Icke, Preston & Balick 1989). 
Model B would therefore also describe 'early' bipolar PNe (e.g. 
NGC 6886, Balick & Frank 2002; see also Hajian et al. 2007) 
which have a barrel-shaped or convex waist against the estab- 
lished convention. In such cases the polar outflows are less de- 
veloped and are yet to create the effect of a pinched waist as seen 
in more evolved nebulae (e.g. NGC 7026, Balick & Frank 2002). 
While these outflows could simply be added to Fig. [2] we have 
omitted them to emphasise that even some more evolved nebu- 
lae will have lost their polar outflows to leave a 'bipolar core'. 
In most cases Model B nebulae are classified possible bipolar 
(PB) since an elliptical model can reproduce the same morphol- 
ogy. If polar outflows are also present then a likely bipolar (LB) 
classification is assigned. 

Multiple-shell (MS) nebulae are also identified in our sample 
(e.g. Chu, Jacoby & Arendt 1987). A separate model is not re- 
quired for these nebulae which are discussed later in terms of 
features drawn from both models (see Sec. |4.3t . Less certain 
cases may be classified as irregular (Irr) typically when there are 
insufficient morphological details available to facilitate a clear 
classification. 

Even though we lack detailed kinematic information for each 
nebula, there is much to be gained from qualitative analysis 
of high-quality images within the framework of these models. 
Apart from the identification of unambiguous morphologies and 
an overall picture of the sample, the models can be applied to 
similar nebulae to unlock a common intrinsic morphology dif- 
ficult to infer otherwise. This approach is also highly comple- 
mentary to kinematic modelling that sometimes require similar 
interpretation to work around limitations inherent in the method 
(e.g. low-inclination and low-velocity situations). Before con- 
ducting our qualitative study we first describe our nomenclature 
regarding low-ionisation structures belonging to Model A. 



3.2. Low-ionisation structures 

At least 17% of PNe exhibit intricate low-ionisation structures 
(LIS) that are often considerably smaller than the overall neb- 
ula diameter (Corradi et al. 1996). Narrow-band [Nil] or [SII] 
images are useful in revealing LIS, but they are best isolated 
by constructing an Hq-+[NII]/[OIII] image that removes density 
variations to highlight just changes in ionisation. Goncalves et al. 
(2001) reviewed PNe exhibiting LIS, their kinematics and their 
purported formation mechanisms (see also Goncalves 2004). 

Studies of LIS have not been helped by a largely inconsis- 
tent approach to their nomenclature in the literature. Perhaps 
the most ambiguous case being the highly collimated 'jet-like' 
pairs often quoted as 'jets' (e.g. Exter et al. 2003) even though 
their velocities not substantially different from their nebulae 
(Goncalves 2004). We describe knots as unresolved or small LIS 
with an aspect ratio of unity, and filaments as knots with an as- 
pect ratio larger than unity (Goncalves et al. 2001). We adopt a 
stricter definition of jets to be opposing polar outflows located 
well outside the main nebula and with typically lower surface 
brightness than equatorial filaments (e.g. Mitchell et al. 2007; 
small red torii in Model A). Jets may appear to be located within 
the radius of the main nebula because of projection effects, with 
their radial position dependent on the inclination of the nebula 
and the height of the jets from the equatorial plane. We avoid the 
'jet-like' nomenclature often applied to highly collimated equa- 
torial filaments that have no relation to our definition of jets. 

4. Nebular Morphologies 

4. 1 . Canonical bipolars 

The canonical bipolars are described by Corradi & Schwarz 
(1995) to have a well-defined equatorial waist or ring-like struc- 
ture from which two fainter, extended lobes diverge in the po- 
lar direction. As a group the majority of them are evolved and 
so are best represented by Model A or some slight variation of 
Model A at intermediate-high inclination. At inclinations lower 
than i ~ 45° kinematic modelling is required to confirm an in- 
trinsic bipolar morphology. When the kinematics exist for low- 
inclination nebulae confirming bipolarity we retain a likely bipo- 
lar morphology (LB) such that our canonical bipolar (CB) clas- 
sification is compatible with the literature. 

Amongst the known sample only NGC 2346 (Bond & Livio 
1990) and A 63 (Mitchell et al. 2007) have been identified as 
canonical bipolar nebulae with inclinations greater than i ~ 45°. 
Such a small fraction is puzzling considering the expectation of a 
high density contrast between the equatorial and polar directions 
for post-CE nebulae (Sandquist et al. 1998). This landscape has 
now changed dramatically with the more representative OGLE 
sample containing at least four canonical bipolars, all of which 
have significantly shorter periods than NGC 2346. These are de- 
scribed below in addition to A 41 which upon reappraisal ex- 
hibits a canonical bipolar morphology. Further bipolars may be 
found in Sec. 14.21 and Sec. 14.41 but many of these require full 
kinematic modelling. 

4.1.1. M 2-19, H 2-29 and A 41 

Figure [3] depicts Gemini GMOS colour-composite images of M 
2-19 and H 2-29 alongside reproduced Hq-+[NH] images of A 
41 (Pollacco & Bell 1997) and PN G126.6+01.3 (Mampaso et 
al. 2006). Apart from NGC 2346, the most canonical bipolar 
morphology of all post-CE nebulae is demonstrated by M 2- 
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Fig. 2. Model A is a simple bipolar adapated from Pollacco & Bell (1997). Model B is a rudimentary adaptation of the Frank et al. 
(1993) model. The inclination angle of the nebulae range from 90° (edge-on) to 0° (pole-on). 



19. It best matches Model A at i ~ 45° and as expected satis- 
fies Gatley's rule (Fig. [TJ. A prominent equatorial ring is seen 
amongst all nebulae in Fig. [3] and this is complemented by a 
separate inner nebula tilted to the major axis in all except H 
2-29. Both A 41 and M 2-19 are the first confirmed post-CE 
nebulae to show this quadrupolar trait consistent with preces- 
sion associated with a binary nucleus (Manchado et al. 1996; 
Mampaso et al. 2006). A binary nucleus appears to be present 
in PN G126.6+01.3, but the definitive proof of an orbital period 
has not yet been realised. In M 2-19 the feature is only visible 
in the unsharp masked [OIII] image as included in the colour- 
composite. Inspection of Model A reveals an inclination effect to 
be likely responsible for the less-developed lobes of the eclips- 
ing H 2-29 compared to M 2-19 and A 41. Even more evolved 
than M 2-19 and A 41 is PN G126.6+01 .3 which exhibits jets ex- 
tending two arcminutes either side of the central star (Mampaso 
et al. 2006). H 2-29 also shows some tentative evidence for a 
polar outflow in the S W direction (Fig. [TJ as well as LIS partic- 
ularly along its edges. 



4.1.2. H 1-33, PHR 1757-2824, PHR 1759-2915 and PHR 
1804-2645 

The remaining canonical bipolars are both highly collimated 
(H 1-33 and PHR 1757-2824) and less collimated (PHR 
1759-2915 and PHR 1804-2645) nebulae (Fig.Q]). Gatley's rule 
applies as expected for H 1-33 and PHR 1804-2645 with band 
4 detections in the waist region of both nebulae. All objects ex- 
hibit eclipses except for H 1-33 whose inclination may be similar 




Fig. 3. Gemini GMOS colour-composite images of M 2-19 and 
H 2-29 made from Ha+[NII] (red), [SII] (green) and [OIII] 
(blue) reveal equatorial rings. Rings are also seen in repro- 
duced Ha+[NII] images of A 41 (Pollacco & Bell 1997) and 
PN G126.6+01.3 (Mampaso et al. 2006). Image dimensions are 
30 x 30 arcsec 2 and some unsharp masking has been applied to 
all except H 2-29. 
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Fig. 4. Gemini GMOS colour-composite 1.7 x 1.7 arcmin 2 im- 
age of Sab 41. Colours are Ha and [Nil] (red), [Nil] (green) and 
[OIII] (blue). 



to M 2-19 at i ~ 45°. A bipolar waist is apparent in both PHR 
1759-2915 and PHR 1804-2645 along PA~ 90° from which 
very faint lobes diverge north and south. This is less evident in 
PHR 1759-2915 which seems to be a fainter version of PHR 
1709-3931 (Parker et al. 2006). 

4.2. Likely bipolars 

4.2.1. Sab 41 

Discovered by Zanin et al. (1997), Sab 41 was first imaged by 
Cappellaro et al. (2001) to reveal an elliptical nebula measur- 
ing 42" x 77". It was rediscovered independently by Parker 
et al. (2006) as PHR 1748-3538 who noticed faint extensions 
beyond the main nebula in the considerably deeper SHS image 
(Fig. [I). Figures |4] and |5] present our GMOS images which pro- 
vide a much more complete view of the nebula and its exten- 
sions. The main elliptical ring measures 30" x 50"in [OIII] and 
slightly larger in [Nil] at 30" x90". The [Nil] ring and in partic- 
ular its terminating S W end is seen amongst bipolar nebulae like 
Lo 18 (Schwarz et al. 1992), while the [OIII] emission traces the 
main nebula body. Inside the elliptical ring are several radially 
distributed low-ionisation filaments with extended tails pointing 
away from the central star. Also visible are polar jets that end 
in shocked [Nil] emission 1.55 arcmin from the central star just 
like those seen in A 63 (Mitchell et al. 2007). In the plane of the 
sky these are bent 30° away from the major axis defined by the 
[OIII] outflow. All these features make Sab 41 an exemplary fit 
to Model A at/ -40°. 

4.2.2. NGC 6337 and Sp 1 

Although NGC 6337 is a relatively new discovery amongst the 
known sample (Hillwig 2004; Hillwig, Bond & Afsar 2006), pre- 
vious interest in its low-ionisation structures gave rise to the first 
kinematic study of the nebula (Corradi et al. 1996, 2000). Figure 




Fig. 5. Gemini GMOS image of Sab 41 constructed by dividing 
[OIII] by Ha+[NII]. Note the opposing faint [NII]-rich end-caps 
1 .55 arcmin from the central star marked in this 4x2 arcmin 2 
image. 



|6]reveals these low-ionisation structures in an [OIII]/[NII] image 
reconstructed from archival data (Corradi et al. 2000). Corradi 
et al. (2000) found only a small difference in the measured ra- 
dial velocity either side of the main ring and interpreted this as 
evidence for a pole-on bipolar nebula with low-ionisation fila- 
ments radially distributed around a projected equatorial density 
enhancement. A low-inclination was vindicated by lightcurve 
analysis which gave a best-fit inclination of i = 9° (Hillwig et 
al. 2006), while a bipolar nebula is supported by Gatley's rule 
being satisfied by a ring of Spitzer IRAC band 4 emission (Fig. 

Having established NGC 6337 as a pole-on bipolar befit- 
ting Model A at i — 10°, explaining the pairs of low-ionisation 
features outside the main nebula body is now straightforward. 
Kinematic studies tell us the motion of pair 'a' is slow and pair 
'b' is fast (Corradi et al. 2000; Garcfa-Dfaz et al. 2009). Both 
studies support a polar outflow interpretation for pair 'b' given 
its high velocity, however pair 'a' is more difficult to accurately 
model given its combination of low-inclination and low-velocity. 
Garcfa-Dfaz et al. (2009) recognised the inability of their model 
to constrain the distance of pair 'a' from the equatorial torus, so 
it is still plausible that pair 'a' represents a slow-moving polar 
outflow. The modelling effort is further complicated by a high 
likelihood that both outflows are bent, not only to each other, but 
also to the polar axis. Nevertheless, the identification of essen- 
tially identical polar outflows in A 63 and Sab 41 (small red torii 
built into Model A) lead us to strongly suggest pair 'a' is indeed 
a similar polar outflow. The strong resemblance between Sab 41 
(Fig.g) and NGC 6337 is discussed in Sec. 1331 
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Fig. 6. The low inclination bipolars NGC 6337 and Sp 1. The 
first column shows the [OIII]/[NII] image of NGC 6337 (Corradi 
et al. 2000) and an unsharp masked SHS/SSS colour-composite 
image of Sp 1 . These are accompanied by their Spitzer IRAC 
colour-composite counterparts in the second column. Image di- 
mensions are 2.5 (NGC 6337) and 2.0 (Sp 1) arcmin a side. 



Bond & Livio (1990) similarly suspected Sp 1 (Fig. [6} to be 
a pole-on bipolar because of its ridged appearance. A kinematic 
study has apparently confirmed this classification (Mitchell, priv. 
communication to Zijlstra 2007) and this is supported again by 
Gatley's rule being satisfied by a ring of Spitzer IRAC band 4 
emission as in NGC 6337 (Fig. [6]). 

4.2.3. A 65 and PPA 1759-2834 

Abell 65 has traditionally been classified as an elliptical PN 
(Bond & Livio 1990; Pollacco & Bell 1997), however this clas- 
sification did not take into account deeper [OIII] images that re- 
veal considerably more outer structure (Walsh & Walton 1996; 
Hua, Dopita & Martinis 1998). Figure [7] reproduces the Hua et 
al. (1998) images as a colour-composite image with the inner re- 
gion shown at lower contrast overlaid on the full nebula at higher 
contrast. Walsh & Walton (1996) suggested the dark lane seen 
along the major axis of the bright elliptical component was actu- 
ally a gap similar to what we see in Model A for i ~ 45° aligned 
along PA~ 45°, leading to A 65 being classified as a bipolar 
waist if only the shallow images are examined (e.g. De Marco 
2009). However, adopting this model does not fit the observed 
outer lobes of A 65, nor does it explain the large [Nil] regions at 
either end of the main body. Furthermore, the central star is not 
placed in the centre of the gap (which itself is rather irregular), 
whereas one may expect from the model a more centred position 
and a sharper, more defined gap. Goldman et al. (2004) discuss 
similar dark patches in NGC 1360 and attribute their origin to 
either lower gas densities or higher extinction. We reject the lat- 
ter possibility based on unpublished Spitzer IRAC photometry 
conducted under program ID 40020 (PI: Fazio) which shows the 
dark patches of NGC 1360 are still present in the mid infra-red. 




Fig. 7. A colour-composite image of A 65 constructed using 
deep [Nil] (red), Ha (green) and [OIII] (blue) images taken from 
Hua etal. (1998). 



This must also be the case for A 65 since Walsh & Walton (1996) 
found no change in HafH/3 across its dark lane. 

Taking a fresh approach a good fit to the outer lobes and 
lemon shaped cross-section seems to be Model A at i ~ 10° and 
PA~ 135°. This interpretation accommodates the [Nil] regions 
as polar outflows which may otherwise be difficult to explain. 
At lower contrast the [Nil] regions seem fragmented and this, 
combined with their conical shape, appears consistent with sim- 
ulations of a ballistic (polar) outflow (e.g. Fig. 6 of Raga et al. 
2007). Higher resolution images and kinematics of the [Nil] re- 
gions are needed to confirm this hypothesis. 

Another case in point for deep, multiple emission line imag- 
ing is PPA 1759-2834. In Paper I the low-resolution SHS image 
shows a round nebula, however the GMOS [OIII] image reveals 
a dramatically different picture (Fig.Q]). A round inner nebula is 
resolved to have internal structure suggestive of a Model B fit at 
i ~ 45°, while two faint extensions are seen in the polar direc- 
tion which are absent in the Ha+[NII] image (just as for A 65). 
These extensions lead us to classify the nebula as likely bipolar 
despite the strong Model B resemblance. A brightened NE rim 
is probably caused by ISM interaction, and while this could be 
an alternative explanation for the extensions (e.g. HFG 1, Sec. 
14.5.21 ). this is unlikely since their position is not coincident with 
the brightened NE rim. 

4.2.4. K 1-2 

All previous studies of K 1-2 have focused solely on the mul- 
tiple opposing pairs of low-ionisation filaments (Corradi et al. 
1999; Exter et al. 2003), rather than the overall nebula structure. 
Figure reproduces the ESO NTT/EMMI [OIII] and [Nil] im- 
ages of K 1-2 from ESO program ID 55.D-0550 (Exter et al. 
2003) which have been enhanced to reveal the outer structures. 
The orientation of the images was chosen for best comparison 
with Model A which is a very good match at the inclination 
of ;=40-50° derived from analysis of the nucleus (De Marco et 
al. 2008). The splayed nature of the multiple opposing filament 
pairs (marked in Fig. [8} seem to trace out the equatorial ring as 
in Model A. Low filament velocities common to both K 1-2 and 
NGC 6337 (Corradi et al. 1999, 2000; Goncalves 2004) further 
support this interpretation. In the central region an apparent in- 
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Fig. 8. Enhanced [OIII] and [Nil] images of K 1 -2 measuring 90 
x 170 arcsec 2 . Low-ionisation filaments are traced by red lines. 



ner gap and elongation in the N-S direction closely mirrors the 
elliptical shape and cross-section of the gap in Model A. We in- 
terpret the ionisation structure of the 'ears' or 'multiple-shells' 
referred to by previous authors as shocked outflows similar to IC 
4634 and M 3-16 (Raga et al. 2008; Sec. 1431] ). 

4.2.5. A 46 

With its eclipsing nucleus, comparison of A 46 with Model A 
should be similar to A 63. Figure [9] presents a previously un- 
published WHT Taurus-2 30 minute [OIII] exposure of A 46 
taken on 18 May 1995. Although of similar depth to the Pollacco 
& Bell (1997) image, the WHT image enables clarification of 
some features identified by them. The inner contours best fit 
Model A if we align the polar axis with PA~55°, however it 
is the outer contours that are of most interest. Pollacco & Bell 
(1997) claimed a 'bow-shock' feature towards the NE, however 
we only find an indentation. We attribute this feature to probable 
ISM interaction where the concave entry (NE) and convex exit 
(SW) points mark the path of the ISM flowing through the neb- 
ula (e.g. Fig. 1 of Dgani & Soker 1998). This is consistent with 
the lower surface brightness SW lobe (Bond & Livio 1990) and 
probably explains the bridge-like feature remarked by Pollacco 
& Bell (1997) at PA-135 as just a contrast difference. 

4.2.6. BE UMa, Bl 3-15 and M PA 1759-3007 

BE UMa appears to be clearly bipolar with an equatorial ring 
and polar outflows (Bond 2000). However, if we adopt Model A 
at the inclination of the system derived from the eclipsing nu- 
cleus (84°), then the apparent inclination of the equatorial ring 
seems too low. This may be explained by its more evolved and 
therefore distended nature and we await a more detailed kine- 
matic study of the nebula. We suspect Bl 3-15 to be a bipolar 
PN complete with probable jets as in A 63 that befit Model A. 
An unresolved nebular core does not allow us to narrow down an 
inclination and while one could argue for an elliptical morphol- 
ogy, the satisfaction of Gatley's rule (Fig . |T} supports our bipolar 
classification. Similarly, for MPA 1759-3007 a central overden- 
sity and probable bipolar lobes are present, suggesting Model A 




Fig. 9. Contoured WHT Taurus-2 [OIII] 2x2 arcmin 2 image of 
A 46. 



is a plausible fit for an inclination angle of i ~ 55°, but higher 
quality images are required before we can classify the nebula as 
a canonical bipolar. 

4.3. Multiple shells 

In a dramatic departure from Bond & Livio (1990), we have 
identified four PNe with multiple-shell morphologies (MSPNe), 
namely PHR 1 804-29 1 3, Pe 1 -9, K 6-34 and Hf 2-2. MSPNe ex- 
hibit considerable morphological (Chu, Jacoby & Arendt 1987) 
and kinematic (Guerrero, Villaver & Manchado 1998) diversity 
which naturally gave rise to multiple proposed formation models 
(Chu et al. 1987; Chu 1989). Bond & Livio (1990) favoured mul- 
tiple mass ejections on the tip of the AGB to explain the absence 
of MSPNe in their sample because the CE phase would interrupt 
such a scenario. We now know that the simplest double-shell PN 
requires only a combination of photo-ionisation and wind inter- 
action (Schonberner et al. 1997; Corradi et al. 2003) and this 
favours the existence of at least some MSPNe with binary nu- 
clei. 

4.3.1. PHR 1804-2913 

PHR 1804-2913 is unique amongst our sample with a round 
double-shell morphology (Fig.[TJ. Measuring 8" across the outer 
shell is more than twice the 3.5" diameter of the inner shell and 
the western side appears to be brightened by interaction with 
the ISM. The eastern edge of the outer shell has an intensity 
30^-5% that of the inner shell, far too bright to be a true halo 
which would be considerably fainter (Corradi et al. 2003). We 
measured a similar intensity ratio for the outer shell of M 2- 
2 from the Ha IAC catalogue image (Manchado et al. 1996a). 
Unlike M 2-2, there is no evidence for LIS in the outer shell of 
PHR 1804-2913. 

4.3.2. Pe 1-9, K 6-34 and Hf 2-2 

The remaining MSPNe Pe 1-9, K 6-34 and Hf 2-2 are depicted 
in Fig. [10] Overall this group is quite homogeneous with an in- 
ner shell dominated by [OIII] emission surrounded by an outer 
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Fig. 10. Multiple-shell nebulae Hf 2-2, Pe 1-9 and K 6-34. The 
PNIC colour-composite image of Hf 2-2 (Schwarz et al. 1992) is 
constructed from Ha (red), [OIII] (blue) and their sum (green). 
Gemini GMOS colour-composite images of Pe 1-9 and K 6-34 
were made using Ha (red), OIIIC (green) and OIII (blue) ex- 
posures. The second K 6-34 panel features an Ha-OIIIC image 
that reveals faint [Nil] emission features outside the main nebula 
(overlaid). Image dimensions are 45 x 45 arcsec 2 (Hf 2-2) and 
30 x 30 arcsec 2 (Pe 1-9 and K 6-34). 



partial shell of low-ionisation knots and small filaments. Unlike 
Hf 2-2 and Pe 1-9, the [OIII] emission in K 6-34 does not extend 
to the outer shell. This suggests a different morphology for K 
6-34, namely a likely bipolar given the satisfaction of Gatley's 
rule and a probable jet revealed by [Nil] emission surrounding 
the nebula similar to NGC 6337. 

The O'Dell, Weiner & Chu (1990) model for NGC 2392 
seems most suitable for the sample. This model is equivalent to 
adding the low-ionisation components of Model A to Model B. 
Hf 2-2 and K 6-34 are probably viewed close to pole-on, while 
reconciling the edge-on Pe 1-9 is a little more difficult since the 
main nebular body is split into two halves. It may be that Pe 1-9 
is just the inner shell of the O'Dell et al. (1990) model viewed 
edge-on without the equatorial LIS. In that case the LIS seen in 
the rim of Pe 1-9 have just migrated from the inner shell, except 
perhaps for one bent filament immediately SW of the CSPN. 
Although their prolate ellipsoidal cores do not make them bipo- 
lar, the MSPNe appear to at least be axisymmetric and they do 
exhibit LIS consistent with the other post-CE nebulae. Many 
other PNe show the same ionisation structure including for ex- 
ample NGC 7662 (Balick 1987, 1998), M 2-2 (Manchado et al. 
1996a) and NGC 2392 (Balick 1987). Whether LIS in this con- 
figuration are a reliable indicator of CE evolution is discussed in 
Sec. 1531 




Fig. 11. An SHS/SSS colour-composite 7x7 arcmin 2 image of 
DS 1 showing two newly identified LIS north of the central star 
(arrowed). 



4.4. Possible bipolars 

4.4.1. DS 1 

In Fig.QT]we identify two northern LIS in DS 1 for the first time. 
The outer most appears more collimated with a trail of smaller 
components while the other is more like a knot. The knot was 
missed by Bond & Livio (1990) due to their lower quality image 
where it appears stellar, while the northern more collimated LIS 
is located outside their image. These LIS are real because they 
appear in both the Ha and B j images. We speculate that the east- 
ern loop of DS 1 represents a polar outflow, placing the major 
axis of the nebula at PA = 90° and the LIS are therefore consis- 
tent with being in the equatorial plane. The overall morphology 
however remains elusive, possibly bipolar, and is probably af- 
fected by density effects (e.g A 65). 

4.4.2. HaTr4, PHR 1756-3342 and PHR 1801-2947 

HaTr 4 (Bond & Livio 1990) appears to be best described by 
Model B at j ~ 45°. The low density contrast required to gen- 
erate Model B requires us to classify HaTr 4 as only a possible 
bipolar. Alternatively, HaTr 4 may resemble H 2-29 but the close 
match with Model B seems more convincing and the original un- 
dersampled lightcurve means that it may or may not be eclipsing 
(De Marco et al. 2008). 

Both PHR 1756-3342 and PHR 1801-2947 exhibit rather 
elliptical nebulae but on closer inspection both have two oppos- 
ing overdensities (Fig. [T}. They may be represented by Model B 
at i ~ 30° and we classify them as possible bipolars for now. 

4.5. Irregulars 
4.5.1. M 3-16 

The intrinsic morphology of M 3-16 is difficult to discern but 
there are a number of similarities that can be drawn with its 
closest analogue, the point-symmetric IC 4634 (Guerrero et al. 
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Fig. 12. Gemini GMOS images of M 3-16 and archival HST im- 
ages of IC 4634 show many similarities in both [OIII] (first col- 
umn) and [MI] (second column). The [Nil] image of M 3-16 
was created in the same way as Fig. [5] Image dimensions are 15 
x 15 arcsec 2 (M 3-16) and 20 x 20 arcsec 2 (IC 4634). 



2008). Figure [12] displays GMOS images of M 3-16 alongside 
archival HST images of IC 4634. The [OIII] emission traces the 
overall nebular morphology and especially the inner most ellip- 
soidal nebula aligned with the major axis at PA~135°. Note the 
inner nebula largely disappears in [Nil] and if a classification 
were to be made elliptical might be the most suitable for both 
nebulae. This is supported by an unusual violet Spitzer false- 
colour (Fig. |T} that mitigates against a bipolar classification us- 
ing Gatley's rule since a number of different scenarios can pro- 
duce this false-colour (Cohen et al. 2007). 

Large outflows at the ends of the major axis are just visi- 
ble in the [OIII] image of each nebula, but they become sharper 
in [Nil]. These features are well-reproduced by models as fast- 
moving photo-ionised neutral clumps (Raga et al. 2008). The 
[Nil] image of M 3-16 also reveals low-ionisation filaments and 
some knots which do not display the same high degree of point 
symmetry as in IC 4634. The eclipsing nucleus of M 3-16 tells 
us the outflows are indeed in the polar direction, which is proba- 
bly also the case for IC 4634. The inclination of IC 4634 appears 
to be low and so in a similar case to NGC 6337 the position of 
the jets is not well constrained by the kinematics and they could 
well be polar outflows (Guerrero et al. 2008). If the inclination 
of IC 4634 were increased to match M 3-16 we might suspect 
its LIS to 'fold in' on themselves to appear similar to those in M 
3-16. 

4.5.2. HFG 1 

In the discovery images (Heckathorn, Fesen & Gull 1982) then 
nebula of HFG 1 apparently resembles the multiple-shell PNe 
IC 1295 (Manchado et al. 1996a) and Wray 17-1 (Corradi et al. 
1999). A more favourable explanation of the brightened SE rim 



is provided by an ISM interaction caused by a high proper mo- 
tion. This was first suggested by Xilouris et al. (1996) and was 
largely confirmed by Boumis et al. (2009) who elucidated a 20' 
long tail of emission consistent with the CSPN proper motion 
via deep imaging. The overall morphology is then indeterminate 
having been substantially modified by ISM interaction. A low- 
ionisation filament was first identified and studied in the discov- 
ery paper, although no kinematic information was given. 

4.5.3. SBS 1150+599A 

De Marco (2009) reproduced an HST F658N image of SBS 
1 150+599 A that shows a round nebula with two collimated fil- 
aments. It is peculiar that spectroscopic studies of the nebula 
(Tovmassian et al. 2001; Jacoby et al. 2002) have not recorded 
any low-ionisation species that could be attributed to these fil- 
aments, especially since the F658N filter is centred on [Nil] 
/16584 with only ~3% transmission at Ha. Attributing the fil- 
aments to Ha emission alone is therefore unlikely and we sus- 
pect further IFU or high spatial resolution longslit observations 
may recover [Nil] emission. The filaments therefore would not 
be representative of the intrinsic morphology which is round or 
unresolved. It is interesting to note the Spitzer IRAC colour- 
composite image has the same violet hue as M 3-16 so we cannot 
classify the nebula as bipolar using Gatley's rule. 

4.5.4. BMP 1800-3408, PPA 1747-3435 and NGC 6026 

Our deep GMOS image of BMP 1800-3408 reveals many fil- 
aments arranged in an asymmetric shell quite reminiscent of a 
supernova remnant (Fig. [TJ. Yet the hot star confirmed by our 
GMOS spectrum has the same radial velocity as the nebula and 
so is certainly the ionising source of the nebula. Further study 
is required to establish whether this is a genuine evolved PN. 
The nebulae of PPA 1747-3435 and NGC 6026 (Schwarz et al. 
1992) are too diffuse to confidently identify any features for mor- 
phological classification. High quality GMOS images of PPA 
1747-3435 were taken but did not improve over the SHS im- 
ages given the diffuse nature of the nebula. 



5. Discussion 

5.1. A penchant for bipolarity 

We are now able to assess the fraction of bipolars amongst our 
post-CE sample (Tab. [TJ. Removal of the as yet unconfirmed 
PHR 1759-2915, MPA 1759-3007 and PHR 1801-2947 leaves 
us with 30 nebulae, of which there are eight canonical bipo- 
lars (27%). Accounting for inclination effects and incorporat- 
ing Gatley's rule brings an additional 9-10 likely bipolars that 
contribute to a plausible bipolar fraction of ~60%. This bipo- 
lar fraction now includes the two essentially confirmed bipolars 
NGC 6337 and Sp 1 that we left as likely bipolar for consis- 
tency with the literature. It may even be as high as 65-70% if 
some of the possible bipolars are included. Around 10% of the 
sample are MSPNe, some of which may turn out to be likely 
bipolars viewed pole-on (e.g. K 6-34 which satisfies Gatley's 
rule). Detailed studies of the irregular nebulae that make up the 
remaining 20% of the sample are needed to resolve their status. 

Our qualitative analysis shows for the first time a penchant 
for bipolarity exists amongst the post-CE nebulae. This is a ma- 
jor step forward from the paradigm where the unusual NGC 
2346 was considered to be the only bipolar (Soker 1997, 1998). 
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It seems highly likely that CE evolution produces bipolar nebu- 
lae more often than single progenitors, but future morphological 
studies of nebulae (in particular those around single CSPN) are 
required to build a control sample for a definitive answer. The 
first point of call for such a study would be those nebulae with 
extensive time-series photometry that rule out a close binary. 
This would include the pulsators and non-detections in Ciardullo 
& Bond (1996), as well as the non close-binary CSPN with 'true' 
identifications from sufficiently sampled fields in Paper I. Even 
if high-quality images were available for the latter sample, there 
may be many members distorting the results which have either 
very late-type main sequence or planet-sized companions below 
the detection limit of current surveys or which may have under- 
gone a previous CE interaction (~Sec. l5.31 >. Another option would 
be to refer to large morphological studies of Galactic PNe as a 
reference point, but since these do not adequately treat inclina- 
tion effects any comparison here would be premature. 

The 30% of our sample with canonical bipolar morphologies 
have considerably shorter orbital periods than NGC 2346. Soker 
(1997, 1998) proposed that NGC 2346 avoided the CE phase 
for as long as possible during its evolution, but clearly the new 
canonical bipolars are free of this restriction. Indeed, the more 
numerous and representative OGLE sample fulfils the greater 
parameter space, whereas NGC 2346 occupies only a narrow 
evolutionary space (Soker 2002). Our results are therefore con- 
sistent with recent predictions, but they still invalidate the first 
critical observations of earlier works (Soker 1997, 1998). 

Another application of the canonical bipolars in the OGLE 
sample concerns the sample selection for Tables 3 and 4 of Soker 
(1997). In general, the former contains bipolars with pronounced 
lobes and waists (e.g. M 2-19), while the latter contains more 
'boxy' nebulae reminiscent of bipolar waists but without pro- 
nounced lobes (e.g. H 2-29). We've highlighted M 2-19 and H 2- 
29 to demonstrate again the influence of inclination on what may 
otherwise be construed solely as evolutionary effects. Consider 
Model A (Fig. [2]), those nebulae viewed close to edge-on (H 2- 
29) will often appear more 'boxy' or compact than those viewed 
at intermediate inclinations (M 2-19) despite them both having 
essentially the same intrinsic morphology. 

5.2. A binary formation scenario for low-ionisation structures 

The origin of LIS is uncertain given the considerable diversity 
in their observed properties (Goncalves et al. 2001), not to men- 
tion the sometimes confusing nomenclature (Sec. 13. 21 ). An accre- 
tion disk established via binary interaction has been proposed to 
power jets (Soker & Livio 1994; Nordhaus & Blackman 2006), 
however binary formation scenarios are yet to be explored for 
other LIS. In what follows we propose a binary scenario in gen- 
eral terms after discussing the occurence of LIS in our morpho- 
logical study. 

Table[T]summarises the LIS identified during the morpholog- 
ical study according to Sec. 13.21 namely whether filaments and 
knots are present in the equatorial plane, or if low-surface bright- 
ness jets were ejected in the polar direction. Our approach to 
nomenclature was motivated by the strong resemblance between 
NGC 6337 and Sab 41, which together bridge the divide between 
the jets of A 63 (Pollacco & Bell 1997; Mitchell et al. 2007) 
and the prominent filaments of some other post-CE nebulae (e.g. 
K 1-2). At least 10 nebulae have equatorial knots or filaments 
(33%), around 20% have jets, while both jets and equatorial LIS 
are present in only 10% of cases. Combining the two types gives 
at least 40% of post-CE which have some kind of LIS, more than 
double the estimated rate of occurrence amongst larger PNe sam- 



ples (Corradi et al. 1996). Together with the emerging trends in 
the sample, e.g. NGC 6337 and Sab 41, the relatively high oc- 
curence tends to suggest a binary origin for LIS. 

We propose that LIS are a natural consequence of a photo- 
ionising wind interacting with clumps of neutral material during 
the PN phase which were already deposited during the CE inter- 
action. Gongalves et al. (2009) require LIS to be mostly a mix of 
dust and H2 in order to resolve a large discrepancy in their ob- 
served densities with model expectations. This material is highly 
concentrated towards the orbital plane with only a small amount 
towards the polar direction (Sandquist et al. 1998). Add a photo- 
ionising wind to this mix and we have the appropriate conditions 
to reproduce the slow-moving equatorial filaments and knots in 
post-CE nebulae (e.g. Raga, Steffen & Gonzalez 2005). Polar 
outflows are better described by fast-moving neutral material in- 
teracting again with a photo-ionising wind (Raga et al. 2007, 
2008). Observationally K 1-2 and M 3-16 fit the model with the 
signature 'reverse bow-shock' ionisation structure in their polar 
outflows (e.g. IC 4634, Raga et al. 2008). The counterparts in 
e.g. A 63 and Sab 41 seem more evolved, which may explain 
the lack of detected [OIII] emission in their outflows. Overall 
this picture is in good agreement with the low-ionisation fila- 
ments in the Helix whose IL cores seem to pre-date the PN phase 
(Matsuura et al. 2009). Much more work is required to develop 
the scenario further, e.g. incorporating the role of spiral shocks 
in the spatial distribution of LIS (Sandquist et al. 1998; Edgar et 
al. 2008). 

5.3. Can low-ionisation structures be used to identify a 
previous common-envelope interaction? 

The sample of post-CE nebulae is still too small to deduce any 
robust trends in nebular morphology. All were discovered via ar- 
duous photometric monitoring campaigns which find only 10- 
20% as close binaries (Paper I). Potentially the elevated oc- 
curence of LIS amongst post-CE nebulae may serve as a useful 
tool in identifying other post-CE nebulae. Pre-selection for LIS 
could prioritise objects for photometric monitoring to quickly 
increase the sample of known binary CSPN. In the event where 
no close binary is found, then the nebula may still have been 
produced by CE evolution if we assume the secondary was de- 
stroyed. This is particularly applicable in the case of jets since 
they may be formed either via a stellar dynamo or via destruction 
of the secondary (Nordhaus & Blackman 2006). 

We are planning photometric surveys of more biased sam- 
ples dominated by LIS, but in the meantime we have surveyed 
the literature as part of a preliminary investigation. Since the 
Goncalves et al. (2001) review, the number of PNe with LIS 
has grown considerably thanks to the PNIC that provides mul- 
tiple emission line archival HST and ground based images (e.g. 
Schwarz et al. 1992; Manchado et al. 1996). Particularly useful 
are the HST F658 [Nil] images used to measure nebular expan- 
sion rates with LIS (e.g. Hajian et al. 2007). We have identified 
~70 PNe mostly from the PNIC for our investigation, of which 
about half have CSPN classifications available. 

Table[3]contains a summary of some of the better studied ob- 
jects in the sample. The third column describes the LIS as jets 
(J), filaments (F) or both (JF), and filaments were chosen over 
jets for some objects where low velocities are known (Gongalves 
2004). For images we refer the reader to the PNIC for most neb- 
ulae and Corradi et al. (1999) for Wray 17-1. The fourth col- 
umn shows the central star classification (Mendez 1991; Werner, 
Heber & Hunger 1991; Kaler, Stanghelli & Shaw 1993; Acker 
& Neiner 2003), while the last column indicates whether photo- 
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Table 3. A select sample of PNe identified to have LIS. 
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Fig. 13. Wolf-Rayet features in the AAOmega spectra of Hf 2-1 
(top panel) and M 3-8 (bottom three panels). 



metric variability has been reported or ruled out (Bond & Grauer 
1987; Bond & Livio 1990; Ciardullo & Bond 1996; Handler 
2003). The variability described here is not from a close com- 
panion, but rather reflects pulsations or variable winds. Ciardullo 
et al. (1999) found no wide companion in NGC 6826, NGC 
6543, NGC 2371 and NGC 5882, and only a possible associ- 
ation in NGC 2392. Hf 2-1 is a new discovery made from our 
deep AAOmega (Sharp et al. 2006) observations described in 
Paper I (Fig. [Oil. 

Our preliminary study shows the majority of the sample 
(~66%) to have emission line nuclei of either Of (~33%), [WR] 
(27%) or PG1159 (6%) type. Around 15% are known or sus- 
pected binaries, while the remaining ~18% are non-emission 
line nuclei (some of which may turn out to be binary). Selecting 
for LIS therefore will give a mix of mostly binary and emission 
line nuclei which will require further observations to separate. 
Almost all the [WR] CSPN in the sample belong to the hot [WO] 
type that have more extreme and chaotic LIS covering their en- 
tire nebulae, presumably due to turbulence from the strong [WR] 
winds disrupting pre-existing LIS (Acker et al. 2002). The LIS 



in post-CE nebulae are much milder in comparison making the 
LIS in [WO] nebulae one of the few (if any) discriminating fea- 
tures between nebulae surrounding [WR] and non-[WR] PNe. 
Representative members of this group include NGC 6751, NGC 
5189, NGC 2452 and Hb 4 (Fig. [Hi. Our discovery of the [WO] 
nucleus of Hf 2-1 suggests HaTr 11 and Pe 1-17 may also turn 
out to be members of this class. 

If we assume LIS are indeed only produced by CE evo- 
lution, then the results of our preliminary study seem to be 
consistent with two speculative binary formation scenarios for 
Hydrogen-deficient central stars (De Marco & Soker 2002; De 
Marco 2008). In short, these scenarios invoke either the destruc- 
tion of a low-mass secondary or a final thermal pulse within the 
CE phase (causing a second CE phase) to enhance mass-loss to 
the point where the central star probably becomes Hydrogen- 
deficient. Such evolutionary paths are therefore highly likely to 
produce similar morphological traits amongst PNe surrounding 
single emission line nuclei. 

Unfortunately no direct evidence for a close binary with a 
[WR] component exists, leaving us dependent on only indirect 
observational evidence. The scenarios provide a lever in that they 
make predictions about the dual-dust chemistry phenomenon 
(Cohen et al. 1999, 2002). The existence of dual-dust chemistry 
amongst only H-deficient CSPN is a critical outcome of the sce- 
narios. Perea-Calderon et al. (2009) puts this link in doubt by 
their apparent discovery of some H-rich central stars with dual- 
dust chemistry. This finding however must be treated with cau- 
tion since at least one object amongst their H-rich sample, M 3-8 
(PN G358.2+04.2), is in fact H-deficient (Fig.Q/j]!. Deeper spec- 
troscopy of the remaining H-rich dual-dust sample is required to 
settle their status and to further test the binary scenarios. It is also 
worth keeping in mind that the binary scenarios are still specu- 
lative and there may be another explanation for dual-dust chem- 
istry that does not involve binary interactions (Perea-Calderon et 
al. 2009). Furthermore, Perea-Calderon et al. (2009) claim the 
binary fraction of 10-20% (Paper I) is too low to support De 
Marco & Soker (2002), however this is disputable given that in 
these scenarios the system would be undetectable because the 
secondary is likely to be destroyed. 

Assuming that the dual-dust predictions are correct, the 
strongest indirect evidence for the scenarios would be the simul- 
taneous presence of strong post-CE morphological traits (jets 
and LIS) around a [WR] central star with dual-dust chemisty. We 
have identified Hb 4 as such an object with both prominent jets 
and low-ionisation structures (Lopez, Steffen & Meaburn 1997; 
Harrington & Borkowski 2000) visible in archival HST images 
(Fig. [Hll. This is complemented by a [W03] CSPN (Acker & 
Neiner 2003) that was recently found to have dual-dust chem- 
istry (Perea-Calderon et al. 2009). The low-ionisation filaments 
are consistent with the [WO] nucleus and the bent jets are quite 
strong evidence some binary interaction has occured (Soker & 
Livio 1994). Independent time-series photometry of Hb 4 is re- 
quired to definitively rule out a current binary as the OGLE-III 
data are most likely affected by nebular contamination (see Paper 
I and in particular Fig. 2 of Paper I). 

6. Conclusion 

A core sample of 30 post-CE PNe was collated from the litera- 
ture (15 PNe) and Paper 1(15 PNe) to conduct a qualitative mor- 
phological study. High quality narrow -band optical imaging was 
sourced mainly from our GMOS programs on Gemini South, 
while some use was made of Spitzer IRAC images to reveal 
H2 emission in nebulae which is a strong indication of bipolar 
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Fig. 14. HST images of Hb 4 reveal jets and outer loops (en- 
hanced Ha background) in addition to many low-ionisation fila- 
ments in the central nebula overlaid as a colour-composite image 
of [Nil] (red) and Ha (green) emission. North is to the left and 
East is down in this 30 x 15 arcsec 2 image. 

morphology (Gatley's rule). Morphological classifications were 
made with the assistance of geometric models adapated from 
Pollacco & Bell (1997) and Frank et al. (1993) which assisted in 
accounting for the influence of inclination effects. Reappraisals 
of a selection of previously studied nebulae around binary CSPN 
yielded some new morphological details and helped frame the 
discoveries from Paper I. Spectroscopy of central stars in the 
OGLE sample will be presented in forthcoming papers in this 
series. 

Our main conclusions are as follows: 

- Nearly 30% of the sample are canonical bipolar nebu- 
lae, breaking the previous paradigm where NGC 2346 was 
the only post-CE canonical bipolar. This more representa- 
tive sub-sample has much shorter periods than NGC 2346 
demonstrating that it is indeed possible to produce a bipolar 
nebula without avoiding the CE phase until the final evolu- 
tionary stages (Soker 1997,1998). 

- A very plausible bipolar fraction of at least 60% (possibly 
as high as 65-70%) is reached once inclination effects and 
Gatley's rule are considered. This is the strongest indica- 
tion yet that the morphologies of post-CE nebulae largely 
satisfy theoretical expectations of a high density contrast es- 
tablished during the CE phase. The link would be strength- 
ened if reliable morphological studies of nebulae around non 
close-binary central stars were available. Such an ambitious 
study requires deep multiple emission line images, extensive 
kinematic modelling and consideration of inclination effects. 

- No round, single-shell nebulae are found in the sample. Four 
multiple-shell nebulae come close, with the roundest being 
PHR 1804-2913, but at least one member of this sample 
may be a pole-on bipolar (K 6-34). The existence of these 
rounder nebulae may support the idea that the CE interaction 
can produce spherical PNe, but this can only be settled by 
dedicated monitoring surveys of round PNe. 

- Low-ionisation structures (LIS) appear to be quite com- 
mon amongst post-CE nebulae (40% occurence c.f. ~17% 
in general PN population). The emerging trends of nebulae 
with very similar morphological traits (e.g. NGC 6337 and 
Sab 41) strongly suggest a binary origin for LIS. LIS seem 
confined to either the orbital plane as radially distributed 
knots or filaments, or to the polar regions as (mostly) low 
surface brightness jets triggered by a dynamo (Nordhaus 
& Blackman 2006). A likely binary formation scenario in- 
cludes the distribution of neutral clumps of dust and H2 dur- 



ing the CE phase into the orbital plane that are then photo- 
ionised by winds during the PN phase. Deep [Nil] imaging 
of the post-CE sample may reveal very low surface bright- 
ness jets to be more common than even found in this work 
(e.g. Vaytet et al. 2009). Evolutionary effects can probably 
explain the absence of LIS in some post-CE nebulae. 
- A preliminary investigation into a larger sample of nebu- 
lae with LIS found two thirds of them to have emission- 
line nuclei. A binary origin may therefore be responsible 
for LIS around emission-line nuclei whereby one or more 
CE phases created the identifiably post-CE morphology and 
dual-dust chemistry (De Marco & Soker 2002; De Marco 
2008). The most extreme LIS appear in nebulae surrounding 
[WO] nuclei distinguishing them apart from from all other 
PNe. Indirect evidence in support of these scenarios is pro- 
vided by Hb 4 which exhibits jets and filamentary LIS, a 
[W03] nucleus and dual-dust chemistry. 
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